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Abstract 



^ Unlike c-plane nitrides, "non-polar" nitrides grown in e.g. the a-plane or m-plane orientation 
encounter anisotropic in-plane strain due to the anisotropy in the lattice and thermal mismatch 

with the substrate or buffer layer. Such anisotropic strain results in a distortion of the wurtzite 
CZ3 unit cell and creates difficulty in accurate determination of lattice parameters and solid phase 

group-Ill content {xsoiid) in ternary alloys. In this paper we show that the lattice distortion is 
^ orthorhombic, and outline a relatively simple procedure for measurement of lattice parameters 

^ of non-polar group Ill-nitrides epilayers from high resolution x-ray diffraction measurements. 

We derive an approximate expression for XsoUd taking into account the anisotropic strain. We 
^ illustrate this using data for a-plane AlGaN, where we measure the lattice parameters and 

estimate the solid phase Al content, and also show that this method is applicable for m-plane 

^ structures as well. 
O 

O 9 Keywords: Al. High resolution X-ray diffraction; A3. Metalorganic vapor phase epitaxy; Bl. 

10 Non-polar; B2. Semiconducting III-V materials. 

^ 11 PACS codes: 81.05.Ea, 78.55. Cr, 81.15.Gh, 61.05.cp, 61.50.Ah 

(N 

CN 12 1. Introduction 

\Q 13 Group Ill-nitrides semiconductors have potential applications in optoelectronics and mi- 

14 croelectronics devices. Nitride semiconductors epilayers grown along the (0001) c-axis of the 

15 wurtzite crystal structure suffer from strong undesirable spontaneous and piezoelectric polar- 

16 ization fields, which give rise to internal electrical fields [1] and impair device performance. In 

17 quantum wells these fields spatially separate the electrons and holes reducing the overlap of 
their wave functions, and causing a reduction of the recombination efficiency and red-shift of 

19 the emission peak in light-emitting devices [2]. A solution to avoid the deleterious polarization- 

20 induced electric field effects is to use group-Ill nitride layers in crystal orientations which have 

21 no polarization field in the growth direction, and hence across the device active region [3]. 

22 Therefore, there is extensive ongoing research towards the growth of "non-polar" (1120) a-plane 

23 and (1100) m-plane (± c-axis) group Ill-nitrides. These non-polar a- or m-plane epilayers are 

24 generally grown on r- or m-plane sapphire substrates respectively. The lattice mismatch and 

25 thermal expansion coefficients of these nitride epilayers with respect to the substrate are differ- 

26 ent along || c and _L c-directions. This gives rise to an anisotropic in-plane strain which distorts 
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27 the basal-plane of the hexagonal unit cell. Similar anisotropic differences in lattice mismatch 

28 and expansion coefficients also exist within the different members of the Ill-Nitride family, hence 

29 there is an anisotropic in-plane strain even when growing for example an a-plane AlGaN epi- 

30 layer on a a-plane GaN buffer layer. The distortion of the basal plane has also been observed 

31 for similar reason in case of c-plane GaN grown on a-plane sapphire substrate 

32 In c-plane oriented films, epilayers are under uniform in-plane strain, which deforms the unit 

33 cell but preserves the hexagonal symmetry of the basal plane. In that case the determination of 

34 lattice parameter and estimation of group III molar fraction in the ternary alloy (for example 

35 Al content of AlGaN films) is relatively straightforward [5,f6]. However, for the non-polar orien- 

36 tations, the anisotropic in-plane strain results a distortion of the wurtzite unit cell which breaks 

37 the hexagonal symmetry. Further, the orientation of the a-plane nitride epilayer on r-plane 

38 sapphire substrate results in a reduced number of available lattice points in the reciprocal space 

39 within the limiting-sphere that are accessible for diffraction measurements. This makes deter- 

40 mination of lattice parameters for non-polar nitrides and consequently the estimation of group 

41 III composition of non-polar nitrides very difficult [3 [9] . A few different procedures for lattice 

42 parameter determination of such structures have been discussed in the literature. Darakchieva 

43 et al. [H] detail a procedure that requires measuring several symmetric and skew-symmetric 

44 planes at multiple azimuth positions in an edge-symmetric geometry. Another approach by 

45 Roder et al. pj uses measurements of interplanar spacings derived from a combination of 9 sym- 

46 metric, asymmetric, and skew-symmetric reflections, which are weighted by their corresponding 

47 fwhm values and used in a fitting routine to match to an orthorhombic structure via an error 

48 minimization routine. In this paper we suggest a slightly different procedure for measuring the 

49 lattice parameters for such distorted systems using high resolution x-ray diffraction (HRXRD). 

50 We first show that the anisotropic strain results in an orthorhombic distortion of the unit cell, 

51 and derive a general expression for the interplanar distance dhki in such structures. Using the 

52 interplanar distances determined from multiple symmetric and skew-symmetric reflections, we 

53 obtain the lattice parameters via a standard least-square error minimization routine that is eas- 

54 ily implemented in standard mathematical software packages using a matrix formulation. We 

55 also derive an approximate expression for the ternary alloy composition of Al^^Gai-^^N epilayers, 

56 where the Al content XsoHd is obtained solving the stress-strain tensor taking into account the 

57 anisotropic strain. The procedure is illustrated using measurements on a-plane AlGaN epilayers 

58 grown on AIN buffer layers on r-plane sapphire substrates. We also show that the procedure is 

59 applicable for m-plane nitrides as well. 

60 2. Experiment 

61 The AlGaN epilayers were grown via metal organic vapour phase epitaxy (MOVPE) in 

62 a 3 X 2" closed-coupled showerhead reactor. Trimethylgallium (TMGa), trimethylaluminium 

63 (TMAl), and ammonia (NH3) are used as precursors and Pd-diffused hydrogen (H2) as carrier 

64 gas. About 0.8/im thick AlGaN epilayers were grown on AIN buffer layer. The details of 

65 the growth procedure can be found in Ref. [10]. The lattice parameters measurement were 

66 carried out by using a Philips X'PERT^*^ high resolution X-ray diffractometer with a symmetric 

67 Ge(220) hybrid monochromator and an asymmetric triple-axis analyzer and PIXcel solid state 

68 detector array. The 26 value of a set of planes (1120), (2110), (1210), (1010), (2130), (2110), 

69 (1120), (1011), (1012), (2131), (2132) were measured to confirm the orthorhombic distortion and 
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70 to estimate the lattice parameters. Absorption measurements were done on backside polished 

71 samples using a Gary 5000 spectrophotometer to estimate the band gap of the epilayers. 

72 3. Measurement of lattice parameters 

73 Fig 1(a) shows a schematic diagram of the relative unit cell orientation of an a-plane III- 

74 nitride epilayer on r-plane sapphire. The in-plane epitaxial relationships between the group 

75 Ill-nitride layer and sapphire are [0001]nitride \\ [llOl] sapphire and [1100]nitride \\ [1120] sapphire- 

76 The thermal expansion co-efficients in the respective directions and the lattice mismatch are 

77 shown in the table 1. It is evident that the lattice and thermal mismatch along [0001] and 

78 [1100] are different which gives rise to the anisotropic strain in the overlayer and thus distorts 

79 the basal plane of the unit cell as shown schematically in Fig l(b-c) (solid line). Further, the 

80 thermal mismatch along c (| Actjoooi] I) and along m (| Aajijoo] I) are larger for GaN/AlN compared 

81 to the GaN/sapphire and AlN/sapphire cases. 




Figure 1: (a) Schematic diagram showing the orientation of wurtzite unit cell of (1120) a-plane oriented nitride 
epilayer on (1101) r-plane sapphire substrate. Anisotropic in-plane strain results orthorhombic distortion and 
requires at least three lattice parameters as (sidewall), at (top) and c for complete description of the unit 
cell, (b-c) The dotted lines and solid lines show the basal plane of a perfect and distorted hexagonal unit cell 
respectively, (b) A choice of coordinate axis that involves four independent lattice variables a, b, c and 7, whereas 
in (c) symmetry allows the reduction of one parameter (a—b) and gives a simpler expression for d^ki and also 
helps to visualize the lattice points in the reciprocal space. 



82 3.1. Choice of coordinate axis 

83 For a perfect hexagonal unit cell the inter-planar distance dhki (or dhku where i = —h — k) 

84 between the (/i/cZ)-planes is given by 1/dlf^i = 4/3. [(/i^ + k'^ + hk)/a'^] + /c^. This expression 

85 cannot be used for the distorted structure shown in Fig.l(b or c) (solid line). To obtain a 

86 relatively simpler expression for new dhku a proper choice of coordinate axis is helpful. In 

87 the first choice [Fig 1(b)], the inter-planar lattice distance involves four independent lattice 

88 variables a, 6, c, and 7, so the expression of the d^ki will be complicated. But for the second 

89 choice [Fig. 1(c)], the symmetry in the distorted basal plane allows one to use a=h and reduce 

90 one parameter thus simplifying the expression for the dhki- Further, such choice of coordinate 

91 axis helps to visualize the lattice points in the reciprocal space, as discussed later in section 3.3. 
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Table 1: Thermal expansion and lattice mismatch 





Thermal expansion 




GaN 


AIN 


Sapphire 








coefficient 


(lo^ir-i) 












<^[0001] 


3.17 


5.27 








tt[1100] 


5.59 


4.15 








"[1101] 






4.7 






"[1120] 






4.5 




Lattice 


% w.r.t. GaN 


% w.r.t. AIN 


% w.r.t. Sapphire 


mismatch 


along c 


along m 


along c 


along m 


along c 


along m 


GaN 








4.09 


2.45 


1.01 


16.06 


AIN 


-3.92 


-2.39 








-2.87 


13.29 



92 For the most general unit cell i.e. a triclinic structure, (Lattice sides a ^ h ^ c and angles 

93 a 7^ /3 7^ 7) the inter-planar lattice distance can be described by [12] 



V = ahc\/ 1 — cos^ (X — cos^ /3 — cos^ 7 — 2 cos a cos /3 cos 7 

O 1,2 2 • 2 

Oil = c sm a 
5*22 = o?(? sin^ /3 
5*33 = 0^1? sin^ 7 
5*12 = afec^ (cos a cos /3 — cos 7) 
'S'23 = a^6c(cos /3 cos 7 — cos a) 
S'13 = a6^c(cos /3 cos 7 — cos a) 



94 For a hexagonal lattice with orthorhombic distortion, we can put a = /3 = 90° and a 

95 and thus obtain a simplified form 

1 /i^ + A;2 - 2/ifccos7 / 



4fcz (a sin 7)2 

96 The above expression for d^^^i is similar to that for an undistorted hexagonal structure, with 

97 additional terms cos(7) and sin(7) which take into accounts the distortion of the basal plane. 

98 3.2. Accuracy in dhkn. estimation 

99 Throughout the experiment, the 1Q value for each plane was measured in the triple axis 

100 geometry and the interplanar distance d^kn was estimated from the Bragg condition given by 

101 (i/jfcj/ = A/2 sin 6'/ifcj; . If is the error in the determination of peak position then the corre- 

102 sponding error in estimation of d is given by Ac? = —d cot 9 /\9. The error can be minimized 

103 by careful optimization of the Eulerian cradle as discussed in [13]. To check the consistency in 

104 our measurements we repeated the measurement of the 29 value for the (1120) refiection several 
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105 times by taking out the sample and reloading it. It was found that the change in 2^-value occurs 

106 in the third decimal place, which would lead to a corresponding change in the d value, Ad/d, 

107 of order 10~^. All our measurement arc based on the accuracy of dhki-valucs up to the 4th 

108 decimal place. Further, we have not incorporated a refractive index correction, and changes in 

109 d/ifer values due to temperature fluctuations during the measurement (within 5° C) which both 
no result in changes of Ad/d, of order 10~^. 

111 3. 3. Confirmation of orthorhombic distortion 

112 We outline below a quick way to confirm the nature of distortion of the unit cell to be 

113 orthorhombic. We choose the symmetric (1120) plane and its equivalent planes (2ll0) and 

114 (1210) which are easily accessible in skew-symmetric geometry by adjusting the rotational axis 

115 (f) and i/j {(f) is the angle of rotation about the normal to the sample mounting surface and ijj 

116 is the tilt of the diffracting plane out of the diffractometer plane). From the symmetry of the 

117 wurtzite basal plane we expect that (iii2o='^2iio='^i2io ioi an undistorted hexagon as shown in 

118 Fig. 2a. Our measurements on a-plane GaN show that (inso ^ c?2iio='^l2io schematically 

119 shown in Fig. 2b(Data Shown in Table. 2). This imphes that the unit cell gets compressed 

120 along m-direction and elongated along a-direction. Hence the angle between the axis vector xi 

121 and X2 is less than 120°. For a-plane epilayers the direct measurement of (iiioo and (iiioo is not 

122 possible because of geometry limitations of our diffractometer. 

123 Further, the two skew-symmetric planes (1010), (0110) give identical 2^- values (i.e. dioio = 

124 doiio)- Since dioio is the interplanar distance between BC and AO planes (Fig. 2), and doiio is 

125 the same between AB and OC -planes, so we can write OA=OC or a = b. We next choose a set 

126 of asymmetric planes which involve the c lattice parameter, like (1122), (1122); (1011), (1011), 

127 (0111), (Olii); (2112), (1212); we found that each set of asymmetric planes gives identical 

128 2^- value (within measurements error: ± 0.001°). Further, it is also found that the peak 

129 positions of the (1010) and (0110) have the azimuthal relationship ^loio = ±180° -|- 0oilO) 

130 similar relationship holds true for (2ll0), (1210) planes also. These measurements confirm 

131 that a = (3 = 90°. From these observations, we conclude that the hexagonal unit cell has 

132 orthorhombic distortion. 



Tal)l( 


> 2; Coiiliruia.tiou of 


ortliorh(juil)ic 


distortion. 


Planes 


Measured value (distorted) 


Undistorted 
(dLu) (A) 


(hkil) 


"^dhkii (degree) 


dhkil (A) 


(1120) 


57.6378 


1.5979 


1.5947 


(2110) 


57.9779 


1.5894 


1.5947 


(1210) 


57.9769 


1.5894 


1.5947 



133 3.4- Least square method 

134 For an orthorhombic distorted hexagonal unit cell, we have derived an expression for the 

135 inter-planar distance dhki (Eq.l). Since this expression has only three independent variables a, 

136 c, and 7, these can be evaluated in principle by measuring the 2^-values of only three reflections of 

137 non-equivalent planes [for example (1120), (1010) and (1011)] from which the lattice parameters 
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Figure 2: The diagram showing the basal plane of (a) perfect hexagonal and (b) orthorhombic distorted wurtzite 
unit cell. Different interplanar distance have been marked by arrows for both the structures, showing the 
inequalities. Reciprocal lattice points of (c) perfect hexagonal and (d) orthorhombic distorted structures. 



138 can be estimated. But to improve the accuracy of measurement it is better to measure multiple 

139 reflections and minimize the error using a least square method. 

140 Here we can write 7 = 120° — 6. If \6\ is small (< 1°), we can write sin(7) = sin(120° — 6) = 

141 1/2(1 - y^6) and cos(7) = cos(120° - 6) = ^3/2(1 - 6/^3). Substituting back to the eq(l), 

142 we can write an approximate expression: 



1 _ h'^ + e + hk{l + P 

4 /i^ + p + hkji + v?,5) / 5 

3 V ~7!; 

4 + P + hk{l + ^5) / 25 \ P 
~ 3 a2 y + c2 

4 (/i2 + fc2 + /ifc) 4 {2h^ + 2^ + 5hk) ^ f 
3 a2 3^ a2 c2 

143 Here we have neglected the higher order (5^-term. By rearranging it, we can write the 
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144 expression in the form of 3 linear variables: 



(P 



{h? + + hk) 



145 
146 
147 
148 
149 
150 
151 
152 
153 



3 

I /a" 



+ 



3v^ 



(2/^2 + 2k'^ + hhk) 



(2) 



Assuming xi = 1/a^, X2 = S/a^ and 2:3 = 1/c^ we can express Eqn.(2) in a linear form 
with variables Xi, X2 and a;3. Now for a set of n reflecting planes say, {hi,ki,li), (/i2,/^2,^2) 
....{hn, kn,ln)] wc wiU gct n-cquations which can be expressed in a matrix form Ap=D, where 
A is a (3 X n) matrix whose clement [A]gj is the co-efficients of the variable Xj {j = 1, 2, 3) of the 
g-th equation {q = 1,2, ...n), p={xi X2 X3), and 'D={l/dl l/dl ... l/rf^J^, where dn = dh„k„ir, 
and T denotes transpose of the matrix. The required matrix, p, whose elements contain the 
lattice parameters can be obtained by solving the matrix equation p=(A^A)~^ (A^D). This 
is of a form that can be easily implemented in standard software packages like Mathematica and 
Matlab. 



3.5. Example 

For our a-plane GaN on AIN buffer layer, we have measured the 2^- value of a set of nine 
reflections as shown in Table 3 and then calculated the matrix A and D. 



Table 3: 2^- value for the set of planes for a-plane GaN 



Planes 


'^dhkii (degree) 


dhkii (A) 


(1120) 


57.64 


1.5980 


(1010) 


32.38 


2.7628 


(2130) 


94.86 


1.0459 


(2110) 


57.98 


1.5894 


(1122) 


69.00 


1.3600 


(1011) 


36.85 


2.4373 


(1012) 


48.13 


1.8889 


(2131) 


97.41 


1.0253 


(2132) 


105.18 


0.9698 



/4.00 


-6.93 


0.00\ 




/0.39\ 


1.33 


-1.54 


0.00 




0.13 


9.33 


-15.39 


0.00 




0.91 


4.00 


-0.00 


0.00 




0.39 


4.00 


-6.93 


4.00 


; and D = 


0.54 


1.33 


-1.54 


1.00 




0.17 


1.33 


-1.54 


4.00 




0.28 


9.33 


-15.39 


1.00 




0.95 


\9.33 


-15.39 


A.OOj 




1^1.06/ 



158 By solving the matrix equation for p using Matlab, we obtain the solution for the lattice 

159 parameters as a = 3.1788 A, c = 5.1774 A, and with 6 = 0.36°, with standard deviations for 

160 cr(a) = 2.9 X 10~^, cr(c) = 6.2 x 10"^ and a{6) = 4.1 x 10"^ respectively. We note that the value 

161 of 6 obtained is small and satisfies the assumption of 5 < 1°. Following the same procedure we 

162 have estimated the lattice parameters for the set of a-plane AlGaN samples, details are shown 

163 in Table 4. 

164 3.6. Procedure for m-plane structure 

165 In case of m-plane nitrides are grown on the m-plane sapphire substrates. Since the lattice 

166 and thermal mismatch along c and a-axis are different, the distortion in the basal plane will be 

167 similar to the a-plane nitrides. Figure 3(a) shows the orientation of unit cell of m-plane nitride 

168 on m-plane sapphire substrate. 




Figure 3: (a) Schematic diagram showing the orientation of wurtzite unit cell of (1100) m-plane oriented nitride 
epilayer on (1100) m-plane sapphire substrate. Anisotropic in-plane strain results in orthorhombic distortion. 
The dotted lines and solid line shows the basal plane of a perfect and distorted hexagonal unit cell respectively, 
(b) A choice of coordinate axis that involves four independent lattice variables a, b, c and 7, whereas in (c) 
symmetry allows the reduction of one parameter {a=b) and gives a simpler expression for dhki- 

169 As in the case of a-plane nitrides, an appropriate choice of coordinate axes can result in 

170 a simpler form of dhki- The first choice (Fig. 3b), the inter-planar lattice distance involves four 

171 independent lattice variables a, b, c, and 7, whereas in the second choice (Fig. 3c), the symmetry 

172 in the distorted basal plane allows one to use a=b and reduces one parameter resulting in a 

173 simpler expression for dhki which is exactly identical to Eq(l). Also the distortion in the unit 

174 cell can be verified using a similar procedure by measuring if the inter-planar spacings follow 

175 diiQQ < c?ioio='^oiio- 

176 This then allows the least square method as mentioned above for a-plane system to be used 

177 to estimate the lattice parameters for m-plane nitrides. 



178 4. Estimation of Al content 

For completely relaxed Al^jGai-ajN films the solid phase Al content {xsoUd) can be estimated 
by measuring either the a or c lattice parameter, subject to the validity of Vegard's law 

ao = x.aA + (1 — x).aG (3a) 
Cq = x.ca + {1-x).cg (3b) 
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179 Here (ac, cg) and {a a-, ca) are the lattice parameters of relaxed GaN and AIN respectively. 

180 Rewriting the above equations 

x= {ac- ao)/(aG - a^) 
x = {cg- co)/(cg - Ca) 

181 Ideally, the x value calculated by using either the a- or c-lattice parameters should be the 

182 same. But in practice for epitaxial AlGaN these values differ. As we have seen in Sec. 3 the 

183 AlGaN epilayers have thermal and lattice mismatches with the buffer /substrate which cause 

184 deformation/distortion in the wurtzite unit cell. For accurate estimation of XsoUd using lat- 

185 tice parameters it is hence necessary to take into account the in-plane strain effect caused by 

186 buffer /substrate. 

187 stress-strain tensor matrix 

188 Considering the X-axis along [1120], F-axis along [1100] and Z-axis along [0001] direction, 

189 the strain-stress relation for hexagonal crystals with a C^^ symmetry can be expressed as 









Cl2 


Cl3 








\ 


1 ^xx\ 


^yy 




Cl2 


Cn 


Cl3 











^yy 


(yzz 




Ci3 


Cl3 


C33 











^zz 















C44 








^xz 


















C44 





^yz 


\<yxyl 




V 














Cm/ 


\^xy/ 



191 
192 
193 



194 
195 



196 
197 



where Cij are the stiffness constants. 



For c-plane nitrides system the crystal is free along [0001]-direction (or Z-axis), so azz = 



and 



C 



13 , 



C: 



~^ ^yy) 



(4a) 



33 



For a-plane nitride epilayers the crystal is free along [1120]-direction (or X-axis), so (Txx = 



and 



Cu 
C 



13 



C 

c 



12 



13 



^yy 



(4b) 



For m-plane nitride epilayers the crystal is free along [1100]-direction (or F-axis), so ayy = 



and 



C 



ezz 



12 



Cu 



^yy 



^ XX ^ 

^Vi '-'13 

198 c-plane AlGaN: effect of deformation on determination of Al content 



(4c) 



We briefly review the well known case of c-plane oriented epilayers [5], where the in-plane 
strain is isotropic which deforms the unit cell, but maintains the hexagonal symmetry. So the 
strain values are 



-yy 



(a - ao)/ao. 



Co)/ Co 
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199 Substituting in eq(4a), we obtain 



(a - ao) , (c - Co) 

h 7 = 

ao Co 

200 Where 7 = 2C13/C33. Substituting the expression for ao and Cq and for small strain values, 

201 Xsoiid can be expressed in terms of lattice parameters as 

^ _ a(cG - c) + 7(aG - a) 
a{aG - ca) + i{cg - ca) 

202 4-3- a-plane AlGaN: effect of deformation on determination of Al content 

203 As discussed in Sec. 3 a-plane nitride epilayers have in-plane anisotropic strain which distorts 

204 the basal plane. Hence, it is necessary to incorporate the effect of anisotropic strain to obtain a 

205 correct expression for Al content. 




Figure 4: The dotted curve showing the basal plane of perfect hexagon and the solid curve same for distorted 
unit cell. Here 1/2 di and ^2 are the interplanar distance of (1120) and (1100) planes respectively. These values 
can be estimated for the expression of d^u in eq.l 



206 Here the strain values are 



^xx = {di — diQ)/diQ 
eyy = {d2 — ^20)/ d2o 
ezz = (c - co)/co 



207 where (ii=2c/ii2o=(2asin7)/ \/2{l — cos 7) and (i2=(iiioo = (a sin 7)/ \/2{l + cos 7) are as in- 

208 dicated in the schematic diagram of the distorted basal plane in Fig. 4. The additional subscript 

209 '0' stands for relaxed AlGaN. Substituting these values into eqn. (4b) we obtain: 
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c — Co di — dio d2 — d2o 
h 7i — ^ \- 72 — ^ = 

Co Ctio 020 

210 Assuming small strain and small distortion {\6\ < 1°), this we obtain 

^ _ (cg - c)did2 + lijdiG - di)cd2 + 72(^^20 - d2)cdi 
(cg - CA)did2 + 7i('^iG - diA)cd2 + 72(C?2G - d2A)cdi 

211 where the additional subscripts 'G' and 'A' stand for GaN and AIN respectively, and 7i=Cii/Ci3 

212 and 72=C'i2/Ci3. For m-plane AlGaN the expression of x has a similar form except with 

213 7i=Ci2/Ci3 and 72=Cii/Ci3. 

214 4-4- Example - solid phase Al composition for a-plane AlGaN 

215 Following the procedure described in the previous section, we have estimated the lattice 

216 parameters and solid phase Al composition for the series of a-plane AlGaN samples grown over 

217 the entire composition range. The details of the lattice parameters derived from HRXRD and 

218 the corresponding Xsoiid calculated are shown in Table 4. This also shows the values of Xgoiid 

219 derived independently from optical transmission measurements 



Table 4: Comparison between x-values obtained from X-ray and transmission. 





a 

(A) 


c 

(A) 


7 

(degree) 


di 

(A) 


d2 

(A) 


solid 

(XRD) 


solid 

(Trans.) 


0.0 


3.1960 


5.1785 


119.64 


3.2132 


2.7628 


0.00 


0.00 


0.2 


3.1783 


5.1507 


119.67 


3.1940 


2.7479 


0.18 


0.19 


0.4 


3.1608 


5.1050 


119.68 


3.1762 


2.7328 


0.38 


0.39 


0.5 


3.1535 


5.0928 


119.70 


3.1654 


2.7269 


0.46 


0.48 


0.7 


3.1362 


5.0532 


119.78 


3.1466 


2.7130 


0.66 


0.66 


1.0 


3.1116 


4.9777 


120.02 


3.1107 


2.6950 


1.00 


1.00 



220 The measurement shows that the value of XsoUd estimated from XRD and optical transmission 

221 agrees within ±2%. The Xsohd is slightly lower than the Xgas because of parasitic reaction between 

222 TMAl and NH3. A detailed discussion of the variation of strain and distortion with Al-content, 

223 and its effect on the microstructure is discussed in references [TU] and [TT], and is not presented 

224 here as the emphasis of this work is to discuss the procedure rather than the results. 

225 5. Conclusion 

226 In conclusion, we have observed that the anisotropic in-plane strain results in an orthorhom- 

227 bic distortion in the wurtzite unit cell for non-polar a-plane nitrides. We have suggested a quick 

228 method for confirming such a distortion to be orthorhombic, and derived an expression for dh^i 

229 value for such distorted unit cells. We have also provide relatively simple procedure for esti- 

230 mation of accurate lattice parameters using multiple reflections and minimizing the error by a 



11 



231 least square method. Since the orthorhombic distortion creates a difficulty for estimating group 

232 III content in ternary alloys, we have presented a technique which estimates the the correct Al 

233 content in a-plane AlGaN films taking into account the effect of anisotropic strain. We have 

234 also shown that this method is equally applicable for m-plane nitrides as well. These procedure 

235 should be valuable to researchers working on a wide range of non-polar Ill-nitride epilayers. 
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